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angle  between  incident  shock  wave  and  reflected  shock  wave  at 
triple  point  (DELTA) 

angle  between  reflected  shock  wave  and  second  Mach  stem  (ETA) 

deflection  angle  of  flow  from  its  original  direction  while 
passing  through  a  shock  wave  (THETA) 

actual  wedge  angle  (THETA  WALL) 

effective  wedge  angle  (THETA  WALL  PRIME) 

angle  between  the  two  reflected  shock  waves  R  and  R'  in  region 
(1)  (XI) 
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second-triple-point-trajectory  angle  (CHI  PRIME) 

angle  between  reflected  shock  wave  and  wedge  surface  (in  RR)  or 
triple-point-trajectory  (in  MR)  (OMEGA  PRIME) 


FROZEN  RR/MR 


f03Q 


o  inoLPQLPoiriom^^iAj'- 
—  ?  jmfnojra-*- 

<? 


FROZEN 


MRCH  NO 


MNIM  01  3Aimi3y  2HOOW 


(QEGJ 


o  in 


o 


(33S/W) 


I  I  I  I  I 


2.50 


O  O  O  LOOOOOOOO 

(n  j  ifl  hou)ooooo 

-  -  -  -^(Mf\IcnariT>f*“0 


1 1 1  \\m\l 


(03S/H) 


133S/H) 


“^0.00  26.00 


tnooooooo 

r^oioooooo 

-<\l<MOO»U>r*  o 


o 

l/> 


EQUILIBRIUM  °  EQUILIBRIUM 


O  u>  Ou»f\JO 

r-  to  <01010  to 


i'oh  ocrse  oo'oe  oo'S2  oo'Q2  oo'st  oo'oi 


IHOdW 


lO  O  COOlOOJ  o 

r»  r-  (D(DimmA 


WWI 


00 '02  I  OO'SOI 


00 '06  00'Si 


00*09 


oo’Sti  oo’oe  oo*si 


OHOdW 


l  1DEG 
85 


DEG  i 


CU2  VIBRATIONAL 
EQUILIBRIUM 


oo’ose  otrosa  oo’Ofis  oo'ooa  00*091 


oo'oai 


Od/£d  QNH  0d/2d 


0?M 


/  / . 


oo ‘ Oh i  oo'oei 

00*001 

00*09 

(03 

oo  o  o 

o  o 

O  IT) 

oo  o  o 

o  in 

o  r* 

^  o  r-  \r>  zr 

m  Co 

(M  — « 

3H  I  fcld  d33WQ 


1  .  _ 

00  09 1  00  OhI  OO'Oei  OO’OOJ 


i030i 


(33S/W) 


FROZEN  RR  5F6  FROZEN 


SHOCK  MRCH  NO 


Od/Ed  GNU  Od/2d 


030) 


(0EG1 


00  8 h  00 


30)  l  di3Hi 


n 


in  o  in  o  m  o 
a*  =r  m  ro  cm  r\j 


i 

| 

| 

1 

i\  ...J 

..  _ _ 

| 

! 

\\  ; 

j 

v\\ 

i 

i  1 

! 

; 

’  ! 

L _  J 

! 

1)  £di3Hl 

- to 

- 20 

- 15 

Otf>9*<OrM~ 

\^\l 

— 

/ 

/ 
. V 

— 

— 1_ 
I 

i 

- J_ 


6  00 '9  00 *e 

(030)  2013H1 


<P\\UI  IV^J 


OO'OM  00 '021  00*001 


00  ’08  00*09  00  *0h  00*02 


(030)  3W I dd  093W0 


in  o  in  o  in  o  m  o 

3 »  <n  m  <\i  <u  ~  — 


00*091 


OO'Ohl  00*021  00*001  00*08  00*09 


00*0h 


00*02 


o  ooooinoooo 
o  ooi/)OM/)3,m(v 


l/>  »<n<M<M' 


FROZEN  MR  SF6  FROZEN 


. V-  - - - 1 - 1 


08  *ti  00  *ti  02*8 


OOHH/EOHU 


u>  oooo 
r-  in  <n<M 


rvi  —  ^  ^ 


00HU/20HU 


2*4. 00  32.00  <10.00  48.00  56.00  °16.00  24.00  32.00  40 

THETfl  WALL  PRIME  (DEG)  THETR  HALL  PRIME 


0.00 


2H3dW 


o 

O 

o 

o 

o 

o 

in 

oooo 

o 

o 

o 

o 

in 

o 

msffxM 

r- 

in 

CO 

<M 

CM 

•m 

o  cd  o  o  inooooo© o 

<m  co  n  in  r-oinooooo 

—  —  —  ~  ^(Mruco^mr*© 


mnim  01  3Aiitnay  shdow 


oooo 

oooo 

=ri/>r-o 


FROZEN  MR  5F6  FROZEN 


O 

o 


WALL  PRIME  (OEG)  THETfi  WALL  PRIME 


AO-A164  §47  TABULAR  AND  OR API CAL  SOLUTIONS  OF  REGULAR  AND  ARCH 
REFLECTIONS  IN  PSEUDO-. .  CU>  TORONTO  UNIV  DOHNSVIEN 
< ONTARIO)  INST  FOR  AEROSPACE  STUDIES  T  C  HU  ET  AL. 
UNCLASSIFIED  JUN  85  UTIAS-282-PT-2  AFOSR-TR-85-1231  F/Q  29/4 


I  BRAT 
OUILI 


V  I BRRT I ONRL  RR/MR 
EQUILIBRIUM 


OQHU/IQHU 


VIBRATIONAL  RR  „  SF6  VIBRATIONAL 


1030) 


oo'oh?  oo’oie  oo’oei  oo’osi  oo’oei 


oo*o 


Od/£d  ONH  Od/ej 


DEGl 


£H3bW 


SF6  VIBRATIONAL  MR  M  '°EG' 


(33S/H)  Sti 


SHOCK 


1 1 BR 
:ou  i 


10EG) 


CH  NO 


ocros 


OO'Ofi 


oo*oe 


oo  #o 


(030) 


£di3Hi 


o 

UJ 

o  in  o  in  o 
^  zt  m  cn 


momom=rrooj  — 
<\jr\j  —  — 


IBR 

OUI 


10.00 


2.50 


.20,  Ms 


ID  O  O  O  OOOO 

D*  O  if)  O  OOOO 

rv  ru  ro  =a>  i/ir^  o 


MNIM  01  3AI1013U  ?H3bW 


00‘t 


93Q)  £bi3Hi 


o  o  o 
o  on  o 


StilBHl 


ARGON  FROZEN 


ARGON  FROZEN 


ARGON  FROZEN’ 


FROZEN  OMR  H I R  FROZEN 


MACH  NO. 


FROZEN 


CH  NO 


FROZEN 


d2bl3Hi 


imo: 

.duaii 


1030) 


BRIUM 


o 

N- 

OJ 

I 

ll 


GAJI/XJjf^OOO* 

i'  \\™\M 

r~. ' — ; — > — mn - r 


OEGt 


-02  VIBRATIONAL  OMR 
EQUILIBRIUM 


Q  OCPCOr-uDLO  fM  o 
w  ina'OPirrrrr  zr  zr 


r- 

o'. 


ir- 

(*> 


dSIHd 


FROZEN  DMR  CO  d  FROZEN 


FROZEN  DMR  IK,J  FROZEN  DMR 


MRCH  NO.  Me  SHOCK  MACH  NO 


DEG  l 


(33S/W) 


(/) 


o 


ONflL  DMR 
RIUH 


VIBRATIONAL  DMR 
EQUILIBRIUM 


dtiU13Hi 


d2d!3Hl 


/  r1 


FROZEN  DMR  SPb  FROZEN  DMR 


(33S/W) 

o 

Ul 

o  cr* 

00 

r>. 

to 

tn 

dSMDdW 


H  (DEG! 
^  M3 


9 


OOVS 


00  -fi5 


oo-  IS 


00  *  8  ti 


(33S/W) 


FROZEN 


cvj  i/>a)r^ooa> 


WWI 


00  ‘  h2 

00  "02 

00*91  00*2! 

3 

U.I 

Q  O)  ©  r>  a?  i 

dht)13Hl 


d2t)13Hl 


310 


VIBRflTIONRL 


dS013Hi 


APPENDIX  A 


ACTUAL  SIDEWALL  PRESSURE  HISTORY  AND  NUMERICAL  RESULTS 


Recently,  Deschambault  (Ref.  II )  measured  the  pressure  history  on  and 
above  a  40°  wedge  for  the  four  types  of  pseudo-stationary  oblique-shock¬ 
wave  reflection  in  air.  The  numerical  results  generated  by  the  present 
work  are  compared  with  the  experimental  pressure  histories  recorded  by 
Deschambault  and  the  results  from  the  numerical  simulation  done  by  Glaz 
(Ref.  12)  for  a  case  of  DMR  in  air.  The  comparison  is  applied  only  to  the 
regions  around  the  two  triple  points  due  to  the  limitations  of  the 
three-shock  theory.  Hence,  only  the  sidewall  pressure  histories  are 
considered.  The  wedge  surface  and  sidewall  pressure  histories  can  be  found 
in  Ref.  11.  Figure  A- 1  illustates  the  sidewall  pressure  gauge  positions 
used  in  the  experiments  by  Deschambault. 


Initial  conditions  :  Ms  =  3.72,  0W  =  40°,  PQ  =  45  torr  (6  kPa)  and 

T q  =  21.94°C 

The  infinite-fringe  interf erogram  is  presented  in  Fig.  A-2.  Since  the 
flow  is  self-similar,  the  position  of  the  gauge  "path"  across  the  wave 
system  can  be  traced  by  superimposing  Fig.  A-l  onto  Fig.  A-2.  The  traced 
positions  of  the  gauges  in  the  reflection  system  are  shown  schematically  in 
Fig.  A-3.  In  this  run,  both  gauges  6  and  7  were  below  the  first 

triple-point  trajectory.  Gauge  6  measured  the  pressure  jumps  across  the 
Mach  stem  M  from  state  (G)  to  state  (3)  close  to  the  wedge  surface,  across 
the  rolled  up  slipstream  and  across  the  tip  of  the  second  Mach  stem  M'  to 
state  (5).  Gauge  7  is  located  slightly  above  gauge  6  and  it  recorded  the 
pressure  across  the  primary  Mach  stem  M,  across  the  upper  portion  of  the 
rolled  up  slipstream  S  and  across  the  second  Mach  stem  M'  from  state  (2)  to 
state  (5).  Gauge  8  was  located  above  the  first  triple-point  trajectory  but 
below  the  second  triple-point  trajectory.  Gauge  8  was  able  to  measure  the 
pressure  jumps  across  the  incident  shock  wave  I  to  state  (1),  across  the 
primary  reflected  shock  wave  R  to  state  (2),  and  across  the  second  triple 
point  T'  to  state  (5).  Gauge  9  was  positioned  above  the  first  and  second 
triple-point  trajectories.  It  recorded  the  pressure  jumps  across  the 

incident  shock  wave  I  to  state  (I)  and  across  the  second  reflected  shock 
wave  R'  to  state  (4)  just  behind  the  second  triple  point  T'.  The 

experimental  sidewall  pressure  histories  are  shown  in  Fig.  A-4.  The 

measured  pressures  of  the  various  flow  states  are  indicated  in  Fig.  A-4  and 
the  various  states  in  Fig.  A-7. 

The  corresponding  pressure  histories  at  the  four  gauge  locations 
given  by  the  numerical  simulation  of  Glaz  (Ref.  12)  are  presented  in 
Fig.  A-5.  Since  each  station  is  merely  a  mesh  point  on  the  numerical  grid 
network,  the  time  resolution  of  the  numerical  simulation  is  much  greater 

than  that  of  the  pressure  gauge.  In  addition,  the  numerical  records  have 
infinite  rise  times,  whereas  physically  the  gauge  has  a  finite  dimension 
and  rise  times  limited  to  about  4  ps.  The  pressure  histories  generated  by 
the  numerical  code  trace  the  local  and  instantaneous  pressure  variations 
ideally  at  a  point  station.  Experimentally,  a  pressure  gauge  provides  the 
average  pressure  over  its  finite  recording  area  and  the  peak  signals  are 


A.  1 


smoothed  out  because  of  the  limited  transient  response  time.  The  pressure 
jump  across  each  shock  wave  is  indicated  in  torr  and  kPa  units  on  the 

plots.  The  lower  pressure  experienced  at  station  6  than  at  station  7 

across  the  rolled  up  slipstream  is  clearly  shown  in  Fig.  A-5  by  the  sharp 
dip  immediately  behind  the  jump  across  the  primary  Mach  stem  M  at  station 
6  here.  The  pressure  jumps  recorded  by  gauge  8  across  the  incident  shock, 
first  reflected  shock  and  the  second  Mach  stem  are  not  as  sharp  and 
definite  as  shown  by  the  simulated  results.  In  the  measured  pressure 

history,  the  rise  from  P2  to  are  continuous.  However  the  simulation 

has  a  pressure  peak  at  P2,  then  the  pressure  dips  slightly  before  rising 

up  to  P5.  Gauge  8  and  9  both  measured  about  the  same  Pj,  but  gauge  9  has 
a  longer  record  of  P^  than  gauge  8  because  gauge  9  went  from  state  (1)  to 
state  (4)  above  T',  whereas  gauge  8  went  from  state  (1)  through  state  (2) 
to  state  (5)  below  T'.  Due  to  the  limited  time  record  of  the  simulated 

pressure  history  given  at  station  9,  the  shape  of  the  pressure  history  and 
the  value  of  P^  cannot  be  read  and  compared  with  the  measured  value. 
However  P^  can  be  read  out  from  the  isobaric  contours  of  the  entire 
reflection  configuration  (Fig.  A-6). 

The  isobaric  contours  of  the  entire  flow  field  (Ref.  12)  and  the 
enlargement  of  the  flow  regions  around  the  two  triple  points  are  presented 
in  Fig.  A-6.  From  Fig.  A-6(b),  both  P2  and  P3  behind  the  first  triple 
point  are  2  * 1 0 6  dynes/cm2  or  200  kPa  (1500  torr),  and  both  P^  and  P5  behind 
the  second  triple  point  are  3*106  dynes/cm2  or  300  kPa  (2250  torr). 
Besides  extracting  the  pressure  information  from  the  isobaric  contour 
plots,  several  reflection  angles  can  also  be  measured  from  this  figure  for 
comparing  with  the  experimental  values.  The  angles  of  interest  are:  (1) 
angle  between  incident  shock  wave  and  reflected  shock  wave  at  the  first 
triple  point  6,  (2)  angle  between  first  reflected  and  second  reflected 

shock  t|,  (3)  angle  between  first  reflected  shock  and  second  Mach  stem  £, 

(4)  first  triple-point  trajectory  angle  x>  and  (5)  second  triple-point 
trajectory  angle  x'»  Refer  to  Fig.  A-7  for  the  definition  of  these 
angles . 

The  computer  program  that  calculates  the  various  thermodynamic  states 
was  modified  so  that  a  unique  thermodynamic  of  state  could  be  specified  for 
each  flow  region.  This  is  seen  to  be  important  for  doing  analysis  in  air 
where  several  flow  regions  may  still  be  frozen  while  the  rest  of  the  flow 
regions  may  be  in  equilibrium.  Concerning  equilibrium  flow  in  air,  there 
may  be  some  circumstances  where  only  the  oxygen  molecules  are  in 
vibrational  equilibrium,  or  in  some  other  cases,  both  oxygen  and  nitrogen 
molecules  are  in  vibrational  equilibrium.  Since  the  incident  shock  Mach 
number  is  low  in  this  case,  no  dissociation  needs  to  be  included  in  the 
analysis.  To  determine  whether  the  flow  is  frozen  or  in  equilibrium, 
relaxation  length  concepts  are  used.  If  the  relaxation  length  of  an 

internal  degree  of  freedom  is  considerably  shorter  than  a  characteristic 
length,  the  state  is  considered  to  be  in  equilibrium.  Otherwise,  it  is 

assumed  to  be  frozen.  The  vibrational  relaxation  length  of  N2  is 
calculated  according  to  the  data  reported  by  White  and  Millikan  (Ref.  13) 
which  shows  that  below  a  temperature  of  3000  K  at  a  pressure  of  one 
atmosphere,  N2  is  excited  more  rapidly  in  air  than  in  pure  N  2.  The 

vibrational  relaxation  length  of  02  in  air  is  calculated  using  the  relation 
given  by  Lutz  and  Kiefer  (Ref.  14)  for  pure  02.  Also  included  is  the 

correction  factor  introduced  by  Vincenti  (Ref.  15).  This  takes  into 


account  the  experimental  findings  of  Blackman  (Ref.  16)  that  0 y~N 2 
collisions  in  air  are  only  40%  as  effective  as  02~02  collisions  in 
transferring  energy  for  vibrational  readjustment  of  U2‘  According  to  this 
analytical  model,  states  (0),  (1)  and  (2)  are  frozen,  whereas  the  flow  in 

states  (3),  (4)  and  (5)  are  in  02  vibrational  equilibrium.  The  vibrational 
relaxation  length  calculated  by  the  present  program  is  38  mm  in  state  (1), 
16  mm  in  state  (2),  3  mm  in  state  (3),  7  mm  in  state  (4)  and  9  mm  in  state 
(5).  Results  of  this  model  are  listed  in  Table  A-l  under  present  analysis 
(A). 

Glaz  used  in  his  simulation  the  real-gas  equation  of  state  for  air  by 
Deschambault  (Ref.  17),  which  was  originally  developed  by  Hansen  (Ref.  18). 
Glaz's  results  indicated  that  all  flow  regions  of  the  reflection  system 
should  be  in  equilibrium  and  a  frozen  gas  assumption  would  not  be  valid. 
Another  run  was  done  using  the  present  program  based  on  the  assumption  that 
all  five  flow  regions  are  in  vibrational  equilibrium  and  the  results  are 
listed  in  Table  A-l  under  present  analysis  (B).  The  pressure  of  each  flow 
region  and  the  reflection  angles  measured  or  calculated  from  experiment, 
numerical  simulation  and  analysis  are  all  tabulated  in  Table  A-l.  Both  the 
numerical  simulation  (Ref.  12)  and  the  analysis  from  the  present  work  give 
extremely  good  results  for  pressures  in  the  flow  regions  (1)  to  (3)  with  a 
maximum  discrepancy  of  2.5%.  The  measured  of  gauge  7  was  used  for 

comparison  since  gauge  6  gave  a  lower  value  because  the  vortical  flow  of 
the  rolled  up  primary  slipstream  went  past  gauge  6.  Behind  the  second 
triple  point  T',  the  pressures  P4  and  P5  given  by  Glaz  are  consistently 
higher  than  the  measured  values  by  13.6%  and  10.5%,  respectively.  However, 
the  numerical  results  from  the  present  study  gave  similar  results  to  the 
experimental  P4  and  P5  values  with  the  maximum  deviation  only  3.6%  for  set 
(A)  and  7.1%  for  set  (B). 

The  pressures  predicted  by  the  present  analysis  (A)  agree  very  well 
with  the  measured  pressures  from  the  experiment.  However,  the  pressures 
predicted  by  the  present  analysis  (B)  agree  better  with  the  numerical 
simulation  results  than  with  the  measurements  from  experiment. 

The  first  triple-point-trajectory  angle  x  was  predicted  very  well  by 
Glaz  and  the  present  work  with  a  maximum  deviation  of  0.5°.  The  second 
triple-point-trajectory  angle  x'  was  very  well  predicted  by  the  present 
analysis  (A).  Numerical  simulation  and  present  analysis  (B)  give  slightly 
smaller  values  of  %'  than  the  measured  experimental  value.  The  present 
analyses  (A)  and  (B)  give  values  of  the  reflection  angles  6,  £  and  t|  to 
within  1.6°  and  3.5°  respectively  that  of  the  experimental  value  for  this 
case.  Glaz  predicted  the  angles  6  and  q  being  4°  and  3.5°  respectively 
too  high,  whereas  the  angle  £  had  a  value  of  9.5°  too  low.  Present 
analysis  (A)  agrees  best  with  experiment  on  physical  quantities.  Present 
analysis  (B)  give  very  reasonable  physical  angles  except  the  value  of  x' 
is  slightly  low.  Numerical  simulation  predicts  triple  point  trajectory 
angles  well  but  not  so  for  the  reflection  angles. 

A  comparison  of  the  values  of  6  shows  how  well  the  numerical  work 
from  the  present  study  can  predict  the  orientation  of  the  reflected  shock 
wave  R  at  the  first  triple  point  T  relative  to  the  incident  shock  wave  I. 
Similarly,  comparisons  of  the  angle  £  and  q  indicates  the  accuracy  of  the 
numerical  work  in  predicting  the  orientation  of  the  second  reflected  shock 


wave  R'  and  the  second  Mach  stem  M*  relative  to  the  primary  reflected 
shock  wave  K.  The  numerical  simulation  by  Glaz  predicts  a  smaller  angle 
f!  and  a  larger  angle  h  than  the  experimental  and  analytical  values,  and  it 
gives  larger  values  of  P4  and  P 5.  This  is  because  decreasing  Z  will  place 
the  second  reflected  shock  R'  more  normal  to  the  incoming  flow, 
subsequently  increasing  the  pressure  behind  it  in  region  (4). 

Similarly,  as  the  second  Mach  stem  M*  becomes  stronger,  the  angle  n  gets 
larger  and  the  incoming  flow  in  region  (2)  is  now  more  normal  to  the 
second  Mach  stem  M',  thus  resulting  in  a  higher  pressure  P5  in  region  (5). 

Present  analysis  (A)  has  given  results  for  the  pressure  at  all  five 
states  (1)  to  (5)  showing  good  agreement  with  experiment.  It  can  predict 
accurately  the  orientations  of  the  shock  waves  around  the  two  triple 
points.  The  good  agreement  indicates  that  the  application  of  different 
equations  of  state  to  each  region  is  justified  based  on  relaxation  length 
concepts.  The  numerical  simulation  by  Glaz  provides  more  detailed 
information  about  the  entire  flow  field  in  contrast  to  the  present  work. 
However,  his  results  show  larger  deviations  in  the  pressures  behind  the 
second  triple  point  and  in  the  reflection  angles  when  compared  to  the 
experimental  results. 


Gauge  6  is  8.56  cm  from  C 
Gauge  7  is  9.65  cm  from  C 
Gauge  8  Is  7.75  cm  from  C 
Gauge  9  is  8.78  cm  from  C 


Fig.  A- 1  Schematic  diagram  of  gauge  positions  on  the 
specially  constructed  steel  window. 
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=  40°,  If,  =  4  5  torr  (6  kPa)  ,  T0  =  21.94 
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14;.  A-5  Numerical  simulation  of  pressure  history  by  Glaz  (Ref.  12) 
for  the  case  of  double-Mach  reflection  in  air.  Ms  =  3.72, 
0,,  =  40°,  P„  =  45  torr  (6  kPa),  1.  =  21.94  °C. 


Numerical  simulation  of  pressure  contours  by  Glaz  (Ref.  12) 
for  the  case  of  double-Mach  reflection  in  air.  Ms  =  3.72, 
6W  =  40°,  PQ  =  45  torr  (6  kPa) ,  T0  =  21.94  °C  and  y  =  1.4. 
Contours  are  shown  from  0.0  to  3.8  x  106  dynes/cm2  with  an 
interval  of  1.0  x  105  and  labels  are  scaled  by  1.0  x  10_1*. 


APPENDIX  B 


COMPUTER- PROGRAM  LISTING  FOR  THE  ANALYTICAL  SOLUTION 
OF  REGULAR  AND  MACH  REFLECTIONS 


The  computer  program  listed  here  follows  the  method  of  calculation 
discussed  in  Section  2.  The  program  is  very  general  and  is  used  for 
argon,  nitrogen,  air,  carbon  dioxide  and  sulfur  hexafluoride.  It  allows 
a  choice  of  frozen,  perfect  or  vibrational  equilibrium  thermodynamic 
states.  It  can  be  easily  modified  to  include  any  other  gas  by  changing 
SUBROUTINE  GAS. 

The  program  will  first  check  whether  the  input  flow  conditions  can 
give  regular  reflection  by  calling  SUBROUTINE  RR.  If  it  fails  to  obtain 
a  solution,  the  program  will  proceed  to  call  SUBROUTINE  SIMMR  to  calculate 
Mach  reflection  case.  If  the  actual  wedge  angle  has  been  chosen  as  the 
parameter,  the  program  will  predict  a  trajectory  angle  for  the  first 
triple  point.  However,  if  the  effective  wedge  angle  has  been  chosen  as 
the  parameter,  the  program  will  iterate  the  wedge  angle  until  the 
calculated  effective  wedge  angle  is  the  same  as  the  input  one.  After 
solving  for  the  first  triple-point  system,  a  check  of  the  flow  Mach  number 
in  region  (2)  with  respect  to  the  kink  is  done.  If  that  is  a  supersonic 
flow,  the  second  triple-point  system  will  be  solved  by  calling  SUBROUTINE 
TRPL2 .  The  subroutines  are  listed  following  the  main  program 
alphabetically.  Input  data  are  explained  in  the  main  program.  The 
program  is  run  on  the  Perkin  Elmer  3250  Computer  System  at  UTIAS. 
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SET  INITIAL  VALUE  FOR  MACHS 
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IF( IC .EQ.l )  RETURN 

CALI.  GASCOPRSO.DTHPO.DDH.DRHOO.DlSt.DAO.DlW.ISTATO)  PHI1P-PHIHPHI0] 

CALL  GAS(DPRS(2>1)1DTMP(2,1),DUM,DRH01,DUM,DUM.DUM.ISTAT1)  MACH1P-MACH1*( 1 .+SIN( PHI01 )*SIN(PHI01 )/ ( SIN(PH1 1P)*SIN( PHI  IP) ) 
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END 
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